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A Linear Boltzmann Transport model within perturbative QCD is developed for the study of 
parton propagation inside the quark-gluon plasma. Both leading partons and thermal recoil partons 
are tracked so that one can also study jet-induced medium excitations. In this study, we implement 
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fragmentation functions and angular distribution of reconstructed jets. 
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I. INTRODUCTION 

When an energetic quark or gluon propagates through 
a dense partonic system such as the quark-gluon plasma 
(QGP), it interacts with the hot medium through mul¬ 
tiple scattering and induced gluon bremsstrahlung, lead¬ 
ing to transverse momentum broadening and energy loss 
of the propagating parton. As a consequence, spec¬ 
tra of large transverse momentum hadrons and jets in 
high-energy heavy-ion collisions should be suppressed if 
a dense partonic matter is formed in the early stage of 
the collisions. This phenomenon is called jet quench¬ 
ing and has been proposed as a powerful tool for the 
study of properties of the QGP in high-energy heavy-ion 
collisions m • It has been observed in experiments at 
both the Relativistic Heavy-ion Collider (RHIC) [4] [5] 
and the Large Hadron Collider (LHC) jB] as evidence 
for the formation of strongly interacting QGP in high- 
energy heavy-ion collisions. Investigators observed not 
only the suppression of single inclusive hadron spectra at 
large transverse momentum CHID] but also back-to-back 
dihadron m and gamma-hadron correlations U 2 mg. 
The jet quenching phenomenon becomes even more dra¬ 
matic in heavy-ion collisions at LHC when one observes 
the onset of dijets and gamma jets with large asymme¬ 
tries in transverse energy [IMS]. These patterns of jet 
quenching can all be understood quantitatively within 
the picture of multiple scattering and parton energy loss 

mm- 

Within the picture of perturbative QCD (pQCD), the 
interaction between an energetic parton and the medium 
is dominated by small-angle scattering and induced gluon 
radiation. The transverse momentum broadening and en¬ 
ergy loss of the propagating parton depend on the trans¬ 
verse momentum transfer from the medium during each 
scattering which in turn encodes the interaction among 
medium partons at the scale as probed by the energetic 


parton. The jet transport parameter defined as the aver¬ 
aged transverse momentum broadening squared per unit 
length of propagation represents a fundamental property 
of the medium [33] _ 0S1 - Recent phenomenological stud¬ 
ies of single-hadron suppression in high-energy heavy-ion 
collisions within collinear factorized pQCD models give 
extracted values of the transport parameter q ss 1.2 ±0.3 
GeV 2 /fm and q ss 1.9 ± 0.7 GeV 2 /fm at an initial time 
To = 0.6 fm/c for a 10 GeV quark at the center of most 
central Au±Au collisions at ^fs = 200 GeV and Pb±Pb 
collisions at ^fs = 2.76 TeV, respectively 06]. These val¬ 
ues are about two orders of magnitude higher than in 
large cold nuclei [47] 08]. 

Given such large values of the jet transport param¬ 
eter, the energy lost to the medium by a propagat¬ 
ing parton is also very large. Dissipation of this large 
amount of lost energy in the medium can lead to jet- 
induced medium excitations such as supersonic waves or 
Mach cones jH ED]. Study of such jet-induced medium 
excitations can shed light on the transport properties 
of the bulk medium which complements that obtained 
from analyses of anisotropic flows due to collective ex¬ 
pansion. Several types of models can be employed to 
study jet-induced medium excitations such as hydrody- 
namical models [51U54] . string theory [551156] and parton 
transport models 153 ■ The Linear Boltzmann Transport 
(LBT) model [55HDTj combines a kinetic description of 
parton propagation and hydrodynamic description of the 
underlying medium evolution. Since it also keeps track 
of thermal recoil partons from each scattering and their 
further propagation in the medium, the LBT model can 
be used to study both jet transport in medium and jet- 
induced medium excitations. 

Jet-induced medium excitations in the form of ther¬ 
mal recoil partons and their further propagation are also 
important for the study of jet quenching through recon¬ 
structed jets and the corresponding jet profiles. Recon¬ 
structed jets defined through a jet-finding algorithm [62] 
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in experiments are collections of collimated showers of 
hadrons within a jet cone. Neglecting nonperturbative 
effects through hadronization, jets reconstructed through 
jet shower partons and final hadrons are equivalent and 
have approximately the same energy. Since interaction 
in pQCD between an energetic parton and the thermal 
medium is dominated by small-angle scatterings, a large 
fraction of thermal recoil partons from jet-medium in¬ 
teraction fall within the jet cone and therefore should 
become a part of the reconstructed jet. These thermal 
recoil partons within the jet cone can contribute to a 
non-negligible fraction of the total jet energy and modify 
both the jet transverse profile (transverse energy distri¬ 
bution inside the jet cone) and fragmentation functions 
(particle distributions in longitudinal momentum inside 
the jet cone) [50] , In addition, a diffusion wake generated 
behind a propagating parton effectively modifies the un¬ 
derlying fluid background and therefore also affects jet 
reconstruction, the final jet energy, and the jet profile. 
These effects can all be studied within the LBT model 
and are important for understanding the medium modi¬ 
fied jet cross section and jet profiles [6111 . 

In this paper, we describe in detail the LBT model 
that has been developed over the last few years and im¬ 
provements to the early versions that have been used for 
the study of jet-induced medium excitations [551 [55], 
transport, and medium modifications [501 HI]. We also 
employ the model to study elastic energy loss through 
multiple scattering, propagation of jet-induced medium 
excitations, and their contributions to parton distribu¬ 
tion and jet profiles of reconstructed jets using the anti- 
kr algorithm in FAST JET [521 . 

Instead of using the small-angle approximation for par- 
ton scatterings in the early versions of the LBT model, 
we have included the complete set of 2 —>• 2 “elastic” 
parton scattering processes with full scattering matrix 
elements including large-angle scatterings, annihilation 
processes, and flavor conversion. These are necessary for 
more accurate description of parton transport in medium. 
Simulations of parton transport according to the Boltz¬ 
mann equation in the LBT model are based on local 
rates of scattering in a medium that evolves according 
to relativistic hydrodynamic equations. We keep track of 
medium recoil partons from each parton-medium scat¬ 
tering and neglect interactions among jet shower partons 
and thermal recoil partons, wherein lies the meaning of 
“linear” in LBT. The model is therefore Lorentz covari¬ 
ant and can avoid difficulties of implementing the test- 
particle method in a full-fledged parton transport model 
[511J . Furthermore, it also allows us to study jet transport 
and jet-induced medium excitations based on event-by¬ 
event simulations with full fluctuations and correlations. 

This article is devoted to a general description and 
test of elastic scatterings in the LBT model and the cor¬ 
responding jet transport and jet-induced medium exci¬ 
tations in a uniform QGP medium. The final medium 
modification of jets and the effect of jet-induced medium 
excitations will depend strongly on the inelastic parton- 


medium scatterings and the dynamical expansion of the 
medium. These will be described in the future when we 
carry out realistic studies of jet transport and medium 
modification in high-energy heavy-ion collisions. 

The remainder of this paper is organized as follows. 
We first give a detailed description of the LBT model in 
Sec. II. We compare numerical and approximate analytic 
results for local parton scattering rates for different chan¬ 
nels, parton energies, and local temperatures. Differen¬ 
tial scattering rates as functions of the final-state mo¬ 
mentum and scattering angle of the medium recoil parton 
from Monte Carlo simulations are also tested against re¬ 
sults calculated with numerical integrations. In Sec. Ill, 
we study the elastic energy loss from multiple scattering 
and in particular the energy and time or length depen¬ 
dence in a static and uniform QGP medium. Time evolu¬ 
tion of the transverse momentum distribution of the lead¬ 
ing parton is also investigated. Propagation of medium 
recoil partons in terms of jet-induced medium excitations 
are investigated in Sec. IV. A supersonic shock wave 
along with a diffusion wake is shown to arise from the 
propagation of the energy and momentum deposited by 
the propagating energetic parton. In Sec. V, we inves¬ 
tigate contributions from medium recoil partons to the 
reconstructed jet energy using the anti-fc^ algorithm in 
FASTJET jet finding package and study the influence of 
medium recoil partons to the jet energy profile and par- 
ton distributions. Finally, a summary and discussions 
are given in Sec. VI. 


II. THE LINEAR BOLTZMANN TRANSPORT 
MODEL 

Considering only 2 —> 2 parton scattering, jet trans¬ 
port and propagation in a QGP medium can be approxi¬ 
mately described by a set of linear Boltzmann equations, 


f d 3 p 2 

[ d 3 p3 

f d 3 pA 

j {2t:) 3 2 E 2 J 

(2tt) 3 2E 3 j 

1 (2tt) 3 2E 4 


Y [fa(pi)fb(j>2) - fc(P3)fd(P4)} \M ab ^ cd \ 2 
b(c,d ) 

x S 2 (s, t, u){2tt) 4 6 4 (p 1 +P 2 -P 3 - Pi), 

(1) 

where the sum is over the flavor of the initial thermal par- 
ton (6) and all possible scattering channels a + b —>• c + d 
for different flavors of the final partons (c,d). \M ab ^ cb \ 2 
is the corresponding matrix element [65] that is aver¬ 
aged (summed) over the initial (final) spin and color as 
a function of standard Mandelstam variables s, f, and u. 
We have neglected the effect of quantum statistics in this 
study. The phase-space distributions for thermal partons 
in the QGP medium /,; (i = b , d) are the Bose-Einstein 
distribution for gluons and Fermi-Dirac distribution for 
quarks and antiquarks with a local temperature T and 
fluid velocity u = (1, v)/y/l — v 2 . For N c = 3 number 
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of colors in QCD, the spin and color degeneracy fac¬ 
tors are g g = 2 • {N 2 — 1) = 16 and g q = 2 • N c = 6 
for gluons and each flavor of quark or antiquark, respec¬ 
tively. The propagating parton before and after each 
scattering is assumed to follow a classical trajectory with 
fi = (2n) 3 8 3 (p- pl)S 3 (x - Xi - plt/Ei ) (i = a,c) as the 
phase-space density. 

Most matrix elements \M ab ^ cb \ 2 for two-parton scat¬ 
tering diverge at small angle, it, t —> 0, for massless par- 
tons. Such divergences disappear when quasi-particle 
modes due to hard-thermal-loop (HTL) corrections are 
taken into account in the calculation of scattering rates 
and transport coefficients [MUIS]. While an exact formu¬ 
lation of the Boltzmann transport within HTL pQCD is 
difficult, we introduce a screening mass for light quarks 
and gluons in matrix elements for two-parton scattering. 
This is the same approach as in Refs. [66j [67]. This is 
equivalent to introducing a Lorentz-invariant regulariza¬ 
tion condition, 

S 2 (s, t, u ) = 0(s > 2p~ D )9( K —s + p 2 D < t < —fijj), (2) 

in the interaction kernel in the Boltzmann transport 
equations [Eq. 0], where po is the Debye screening 
mass, 


Ed = 




(3) 


for Nf = 3 number of active quark flavors in the QGP. 

Matrix elements for parton scattering we use in this 
paper are leading-order pQCD results. The strong cou¬ 
pling constant a s = g 2 /4-7T will be kept at a constant 
value and should be determined through comparisons to 
the experimental data. Since we will only consider par- 
ton propagation in a uniform medium in this study, we 
set a s = 0.3 throughout the paper. 


A. Scattering rates 


For Monte Carlo simulations of the Boltzmann trans¬ 
port equation, we need to evaluate the rate for a hard 
parton of type a scattering with a thermal parton of type 
b via a specific channel a + b —> c + d, 


1 f d 3 p 2 j" d 3 p 3 I" d 3 p 4 

2E x J (2 tt ) 3 2 eJ (2tt) 3 2E 3 J (2n) 3 2E 4 

X fb{P 2 ) \M\l b ^ cd (s,t,u) S 2 (s,t,u) 
x (2tt) 4 5 (4) (pi +P2-P3 ~Pi)- 


(4) 


Summing over all types of initial thermal partons ( b ) and 
all possible channels with different types of final-state 
partons (c,d), we can obtain the total scattering rate for 
an energetic parton a with thermal partons in a thermal 
QGP medium: 


To = ^ Tabbed. (5) 

b,(cd ) 


These scattering rates for different parton types and dif¬ 
ferent channels will be used in Monte Carlo simulations 
of the LBT model. Shown in Fig. [T] are the total scatter¬ 
ing rates for an energetic gluon (upper panel) and quark 
(lower panel) as a function of the incident energy E in a 
thermal QGP medium with zero baryon chemical poten¬ 
tial at a constant temperature T = 200 MeV. Also shown 
are contributions from different channels. Dominant con¬ 
tributions to the total scattering rate of a gluon come 
from gluon-gluon and gluon-quark scattering. Contribu¬ 
tions from gluon-gluon to quark-antiquark conversion are 
about two orders of magnitude smaller. The dominant 
channel to the quark scattering rate is from quark-gluon 
Compton scattering followed by quark-quark scattering. 
Contributions from quark-antiquark annihilation are also 
many orders of magnitude smaller. 

As shown in Fig. |T[ scattering rates for s-channel an¬ 
nihilation processes decrease with the incident energy as 
T = ap ~ a 2 T 2 /E. Total parton scattering rates are 
dominated by f-channel gluon-gluon, gluon-quark, and 
quark-quark scattering processes, whose rates are ap¬ 
proximately independent of the incident energy E. 

In the high-energy limit E T, f-channel gluon and 
quark scattering cross sections can be approximated by 
their small-angle limits , 


dcr ab 


C ab 


dq 2 ± 

Scattering rates f 
medium are then 


2lTa s ( r - 9 r 

^4 y-'ag ~ ^ 


r 9 « E T 9^9b « 42C A C(3) 

b=g,qi,qi 

r q « ]T v qb ^ qb « A2C F c(3) 

b=g,qi,qi 


qq gj- 

(6) 

i thermal 

QGP 

a 2 T 3 

7 r^D 

(7) 

n 2 T 3 

np z D 

(8) 


where £(3) ~ 1.202 is the Apery’s constant, Ca = N c 
and Cf = {N 2 — 1)/2A C . With the Debye screening 
mass given in Eq. ([3]), the above scattering rates are in¬ 
dependent of the incident energy and proportional to the 
temperature T. As shown in Fig. [l] as solid lines, the 
above analytic estimates of total parton scattering rates 
agree with numerical results very well for E T. We 
also show in Fig. [2] both numerical results (symbols) and 
analytic estimates (solid lines) of total scattering rates 
and contributions from different channels as functions 
of the temperature in a thermal QGP medium for a 
fixed parton incident energy E = 100 GeV. As given 
in the above analytic estimates, total rates as well as 
rates for t channels increase linearly with the tempera¬ 
ture T, while rates for s channels increase quadratically 
with T. Note that the parton scattering rates given here 
are only valid in a medium above the QCD phase tran¬ 
sition temperature T > T c . For a rapid crossover as 
indicated by recent lattice QCD simulations [M], the ef¬ 
fective degrees of freedom around the phase transition 
temperature T c « 155 MeV become ambiguous. One has 
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FIG. 1: (Color online) Flavor-summed scattering rates of dif¬ 
ferent processes r a i)->cd and total scattering rates F a of (a) a 
gluon or (b) a quark (antiquark) as a function of the parton 
energy Eq in a static and uniform QGP medium. Analytic 
estimates for the total rates discussed in the text are shown 
as solid lines. 


to rely on model assumptions to treat parton medium 
interaction H3 linn ed]- 


B. Differential rates 

To sample the initial thermal and final parton four- 
momenta in a 2 2 scattering, we choose the 2 axis 

as the direction of the incident parton a with its four- 
momentum pi = (Ei, 0_L,Pi z )- One can integrate out the 



T (GeV) 



FIG. 2: (Color online) Flavor-summed scattering rates of dif¬ 
ferent processes r a f)->cd and the total scattering rates T a of 
(a) a gluon or (b) a quark (antiquark) as a function of the 
temperature T in a static and uniform QGP medium. An¬ 
alytic estimates for the total rates discussed in the text are 
shown as solid lines. 


four-dimensional 5 function in Eq. Q and express the the 
scattering rate for a given channel in a differential form: 

r ab—»cd (2?r) 4 h J d 0 3 Jd(j) 23 J dE 3 

x f h (E 2 , T ) \M\ 2 ab ^ cd (s, t, u)S 2 (s, t, u) (9) 
E 2 E 3 sin d 2 sin d 3 

Ei(l - cos 6*12) - E 3 (l - cos #23) ’ 
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FIG. 3: (Color online) Energy distributions of the outgoing 
parton in (a) g+b —> g + b(b = g , q, q) or (b) g+g —► q + q scat¬ 
tering processes with different initial parton energy Eq in a 
thermal QGP medium at a temperature T = 200 MeV. Solid 
lines are results from Eq. <§ with direct integration while his¬ 
tograms are from the corresponding Monte Carlo simulations. 


where 

E _ _ £i-E 3 (1 - cos 0 3 ) _ 

2 Ei(l — COS0i 2 ) — £ 3(1 — COS023)’ 

<j>i is the azimuth angle, 9i the polar angle of a parton’s 
momentum pi : and <pij and 0 ^ are the azimuth and polar 
angles between two partons’ momenta pi and pj, respec¬ 
tively. 

One can further carry out the integration over the 
direction of the initial thermal parton b’s momentum 
(^ 21 ^ 23 ) and obtain the final parton c’s angular ( 0 3 ) and 
energy (£ 3 ) distributions. Shown in Figs. [3] and [4] as 



E (GeV) 



E (GeV) 

FIG. 4: (Color online) The same as Fig. [3]in (a) q + q —> g + g 
or (b) q + b —> q + b(b = g, qi, qf) scattering processes. 


solid lines are the energy distributions of an outgoing 
parton in some selected channels of parton-medium in¬ 
teraction in a thermal QGP with a temperature T = 200 
MeV. The gluon-mediunr interaction is dominated by t 
and w-channel gluon-gluon and gluon-quark(antiquark) 
scattering. The outgoing gluon can either be the de¬ 
flected incoming gluon which carries most of the incident 
energy Eq minus elastic energy loss ~ p 2 D /T or the ther¬ 
mal recoil gluon whose average energy is its initial ther¬ 
mal energy ~ T plus the elastic energy transfer from the 
scattering. The energy spectrum of the outgoing gluon 
therefore has two peaks at £ ~ T and Eq as shown in the 
upper panel of Fig. [3] Large-angle and s-channel scat¬ 
terings contribute to the valley region between these two 
peaks which is suppressed by a factor of p 2 D /E 0 T. The 
energy spectrum of the outgoing gluon does not vanish 
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at the incident energy Eq. It instead has an exponential 
fall-off with a slope as given by the temperature, indi¬ 
cating that the energetic incident parton can also gain 
energy due to its interaction with thermal partons. 

The energy spectra of the outgoing quark in the quark- 
antiquark pair production g + g —► q + q (lower panel of 
Fig. i and the outgoing gluon in the quark-antiquark 
annihilation q + q —> g + g process (upper panel of Fig. [4]) 
have a similar structure with two peaks at £ ~ T and E 0 . 
This structure, however, arises mainly from the u and t 
channels of the pair production or annihilation process in 
which the two outgoing partons have an equal probability 
to carry most of the incident energy. Contributions from 
the s channels in these processes are suppressed by a 
factor of /. l 2 d /TEq with a flat energy distribution between 
E~T-E 0 . 

To illustrate the typical energy spectra of thermal re¬ 
coil partons from parton-medium interaction we show in 
the lower panel of Fig. [4] the energy distribution of the 
recoil parton b in q + b —> q + b (b = g, q', q') processes. 
For an energetic quark q , these processes are dominated 
by t channels in which recoil partons carry an energy 
Eb ~ Qt/T- The recoil parton spectra should have a 
power law behavior dN/dEb ~ 1 /q\. ~ 1 /E 2 with a peak 
at E b ~ T + g? D /T. The recoil gluon spectra in the 
s and u channels in the quark-gluon Compton scattering 
q+g —> q+g, however, still has a peak at Eb ~ Eq— fi 2 D /T, 
but is suppressed by a factor of p 2 D /E 0 T relative to the 
t channel. 


Because of the dominance of t and u channels of the 
small-angle scattering, only one of the outgoing partons 
in parton-medium scattering carries most of the incident 
energy for most of the time, whose flavor can be, however, 
different from the incident parton. Therefore, we will 
refer to the more energetic parton in the final state of 
each parton-medium scattering as the leading parton and 
the other less energetic one as the thermal recoil parton. 

For Monte Carlo simulations in the LBT model, we 
change the variable of the energy integration in Eq. Q 
from the outgoing parton c’s energy £ 3 to the initial par- 
ton 6 ’s thermal energy £ 2 . The integrand then can be 
used to sample the thermal parton &’s energy E 2 , polar 
angle 62 and the outgoing parton c’s direction (@ 3 ,^ 23 ) 
with the rejection method. The energy of the outgoing 
parton c £3 and parton d’s four-momentum P 4 can be de¬ 
termined by energy-momentum conservation. Shown as 
histograms in Figs. [3] and [4] are the energy distributions of 
the outgoing parton in some elected channels of parton- 
medium scattering. The agreement between Monte Carlo 
results and those from direct integration in Eq. (|9| pro¬ 
vides numerical verifications of the Monte Carlo code for 
parton-medium scatterings in the LBT model. 

For the purpose of discussing jet-induced medium ex¬ 
citation later, we show in Figs. [5] and [ 6 ] distributions in 
the polar angle of the thermal recoil parton from parton- 
medium scatterings via the same selected channels as 
in Figs. [3] and [4] from both Monte Carlo simulations 
(histograms) and numerical integration of Eq. (solid 
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FIG. 5: (Color online) Polar angle distributions of the ther¬ 
mal recoil parton of different energies relative to the incoming 
parton direction after (a) g + b —¥ g + b (b = g,q, q) or (b) 
g + g —> q + q scattering processes with initial parton energy 
Eq = 20 GeV in a uniform medium at a temperature T = 200 
MeV. Solid lines are for results from the numerical integration 
and histograms for LBT Monte Carlo simulations. 


lines). The polar angle of the thermal recoil parton in 
two different energy ranges, E =0-1 and 1-2 GeV, is 
calculated with respect to the direction of the incident 
parton with Eq = 20 GeV. As shown in the figures, ther¬ 
mal recoil partons with energy significantly larger than 
the typical thermal energy gain their energy and momen¬ 
tum through elastic scattering and therefore are dragged 
along in the same direction of the energetic incident par- 
ton. Soft thermal recoil partons with £=0-1 GeV, on 
the other hand, have large angles, almost perpendicular 
to the incident parton on the average. These features of 
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FIG. 6: (Color online) The same as Fig.[5]for (a) q + q —> g + g 
or (b) q + b —> q + b(b — g, q 1 ,q') scattering processes. 


parton-medium scattering will determine the jet-induced 
medium excitation due to multiple parton scattering as 
we will describe later. 


C. Multiple scatterings 

In the LBT model, we initiate the Monte Carlo pro¬ 
gram by calculating parton scattering rates according 
to Eqs. Q and ([5]) for a range of incident energies 
and medium temperatures which are stored in two- 
dimensional tables for later use. We then approximate 
the propagation of partons in medium with a discretized 
time internal At in the frame of hydrodynamic evolution 
of the QGP medium. We sample the probability for n 
number of parton-medium scatterings within a time in¬ 


terval At according to a Poisson distribution, 

P a , n (E,At ) = [ AJV °(‘ E )]" e —AJT B (g,At) ( 1Q ) 

n\ 

T) ‘ 'll 

AN a (E,At) = At^— T a (p-u,T), (11) 

hj 

where the scattering rate T a (p-u, T) for a parton a is eval¬ 
uated at the local fluid comoving frame with incident en¬ 
ergy E u = p-u and u is the local fluid four-velocity. After 
determining the number of parton scatterings n, one can 
assume these n number of scatterings occur sequentially 
within the time interval At and energy and momentum 
are strictly conserved in each scattering along the classi¬ 
cal trajectory. 

If we choose the value of time interval At to be 
much smaller than the mean-free-path length At -C 
1 /max(r g ,r 9 ), one can approximate the probability for 
at least one parton-medium scattering, 

P a {E,At ) = 1 - e - AAr “( E ’ A ‘), (12) 

as the probability for one parton-medium scattering dur¬ 
ing the time interval At. 

For each parton-medium scattering, we use fractional 
rates T a ^ c d/T a to decide the channel and flavors of 
the initial thermal parton b and outgoing partons c 
and d. The kinematics of the parton-parton scattering 
a + b —> c + dis sampled according to the differential 
rate in Eq. <©• The more energetic parton of the two 
outgoing partons is chosen as the leading parton for the 
next parton-medium scattering. In the LBT model, all 
outgoing partons, both the leading and thermal recoil 
partons, in each parton-medium scattering are recorded 
and are allowed to go through further parton-medium 
scattering in the subsequent time intervals. Each parton- 
medium scattering could in principle accompanied by in¬ 
duced gluon bremsstrahlung. We focus mainly on elastic 
scattering processes in the LBT model in this paper and 
leave the discussion on the implementation of induced 
gluon bremsstrahlung for the subsequent publication. 

To take into account of the back-reaction in the Boltz¬ 
mann transport equation ([!]), we also record the initial 
thermal parton b and its four-momentum p 2 in each scat¬ 
tering process, which we denote as “negative” partons 
and are allowed to transport further according to the 
Boltzmann equation. These “negative” partons should 
be subtracted from the final parton spectra and energy- 
momentum density of the jet-induced medium excita¬ 
tion. Thermal recoil partons and the “negative” partons 
are collectively called jet-induced medium partons in this 
study. 


III. TRANSVERSE MOMENTUM 
BROADENING AND ELASTIC ENERGY LOSS 

In the LBT model, all channels of 2 —► 2 scattering are 
considered for parton-medium interaction in which the 










flavor of the most energetic parton in the final state can 
be different from the incident parton. In the LBT Monte 
Carlo simulations, we designate the more energetic one 
of the two final partons as the leading parton. One can 
then follow the propagation of the leading parton and 
study its transverse momentum broadening and elastic 
energy loss. 


A. pr broadening 


During the propagation of an energetic parton in 
a QGP medium, each parton-medium scattering con¬ 
tributes to the transverse momentum (perpendicular to 
the initial direction of the propagating parton) of the fi¬ 
nal leading parton according to matrix elements of the 
2 —► 2 parton-parton scattering, leading to the increase 
of the averaged transverse momentum squared over time 
or the transverse momentum broadening. The averaged 
transverse momentum squared per unit length is often 
defined as the jet transport parameter q [31]"05, which 
characterizes both the local gluon number density and 
the strength of jet-medium interaction. 

Employing the small-angle approximation of 2 —> 2 
cross sections in Eq. Q, one can estimate the trans¬ 
verse momentum broadening squared per mean-free-path 
length as, 


Qa = (<?i/A)a = 


ps/4 


b,(cd) 


2 ddab—^cd 9 

dq A 


dq 


Pb 


= r a (ql) 


rr 42C(3) 2^3, , a* 

Ca - a s T In — 

7r 4 /j, 


), 


(13) 


where (q^_) is the transverse momentum transfer per 
scattering which is independent of the parton’s flavor, 
C a = Cf = 4/3 for a quark and C a = Ca = 3 for a gluon. 
The variable in the logarithm is defined as In s* = (In s) 
averaged over the kinematics of a single parton-medium 
scattering and s is the center-of-mass energy squared for 
each parton-parton scattering. One can assume, how¬ 
ever, s* = 2cE 0 T with the constant c determined from 
the numerical calculations. 

During multiple scatterings in the LBT model, one can 
follow the propagation of the leading parton and study 
its final transverse momentum distribution and accumu¬ 
lated transverse momentum broadening. Shown in Fig. [7] 
are the transverse momentum broadening per unit length 
from single and multiple scattering of an energetic parton 
with different initial energy E 0 in a static and uniform 
QGP at different temperatures over a length of L = 8 
fm. 

We also show in Fig. [8] the transverse momentum dis¬ 
tributions (solid lines) of the leading parton at different 
times for an initial gluon (upper panel) or quark (lower 
panel) with energy Eq = 100 GeV propagating in a uni¬ 
form medium at a temperature T = 400 MeV. The trans¬ 
verse momentum broadening manifests in these distribu¬ 
tions as the increase of the widths of distributions with 



FIG. 7: (Color online) The jet transport parameter q or trans¬ 
verse momentum broadening squared per unit length of the 
leading parton for (a) an initial gluon or (b) quark with dif¬ 
ferent initial energies going through a single (open symbols) 
or multiple scatterings (solid symbols) in a uniform and static 
QGP medium at different temperatures. Solid lines are an¬ 
alytic results from a single scattering with the small-angle 
approximation. 


time. Such an increase of the width gives rise to an ap¬ 
proximately time-independent jet transport parameter q 
as shown in Fig. [7] To illustrate the evolution of the 
transverse momentum distribution with time beyond the 
change of the width, we also compare each distribution 
with a Gaussian, 


= JL exp(-A) 

dp 2 T (p 2 t ) Pl <p§y’ 


( 14 ) 
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FIG. 8: (Color online) Transverse momentum distributions 
(solid lines) of the leading jet parton at different times for an 
initial (a) gluon or (b) quark propagating in a uniform and 
static medium at a temperature T = 400 MeV. Gaussian fits 
are shown in dashed lines. 


where ( p is the average transverse momentum squared 
from each corresponding transverse momentum distribu¬ 
tion from LBT simulations. The pt distribution from 
a single scattering follows a power law ~ 1 /p^, at small 
angles p^ -C s*/4. At later times, the distribution devi¬ 
ates from the power law behavior due to multiple scatter¬ 
ing. According to Eqs. ([7]) and the total scattering 
rate for a gluon (quark) is T g (T g ) « 1.53 (0.68) 1/fm. 
The leading parton should experience several scatterings 
with medium partons after a few fm/c of time. Its trans¬ 
verse momentum distribution will start to deviate from 


the power law and approach to a Gaussian according to 
the central limit theorem because the transverse momen¬ 
tum transfers during each scattering are independent of 
each other. Since a gluon’s scattering rate is more than 
twice that of a quark, its approach a Gaussian form of 
the pt distribution is also twice sooner than a quark. In 
the high pt region, the power-law-like tail, however, will 
always remain due to a single hard or large-angle scat¬ 
tering. This evolution of the transverse momentum dis¬ 
tribution of the leading parton is also shown in a recent 
perturbative calculation within a HTL pQCD approach 

EH- 


B. Elastic energy loss 


Following the propagation of the leading parton, one 
can also calculate the elastic energy loss in the LBT 
model. Shown in Fig. [9] is the calculated elastic energy 
loss per unit length from LBT simulations as a function 
of the initial parton’s energy for a gluon (upper panel) 
or a quark (lower panel) propagating in a uniform QGP 
medium with a length L = 8 fm at different constant 
temperatures. Also shown in the figure are the simu¬ 
lated elastic energy loss per unit mean-free-path length 
from a single scattering, 


dEl 

dX 


(I», 


(15) 


where v is the energy transfer of the leading parton to 
a thermal medium parton which should depend on the 
energy u> of the thermal parton and the transverse mo¬ 
mentum transfer squared q\ in each scattering. Results 
from a single scattering are in good agreement with that 
from multiple scatterings over a long distance, especially 
for more energetic quarks at lower temperatures. The dif¬ 
ference between results from a single and multiple scat¬ 
terings becomes bigger for lower energy gluons at higher 
temperature, which can be understood by the energy de¬ 
pendence of the energy loss in a single scattering. Using 
the small-angle approximation for elastic scattering cross 
sections in Eq. (|6]) and v k q 2 ^/ 2w, the elastic energy loss 
per mean free path can be expressed as the product of 
the averaged transverse momentum transfer squared per 
scattering {q\) « p 2 D ln(s*/4/i|,) and the 1/w-weighted 
scattering rate [55) : 


/T a 


\ 2uj 


W L )=Ca 


‘X'ra 2 

2 p 


rj-iZ 


D 


(16) 


One obtains the analytic form of the elastic energy loss 
per unit mean free path, 


dEl 

dX 




(17) 


under the small-angle scattering approximation. This 
analytic result agrees with the full HTL pQCD result 
El El in the leading logarithmic approximation and 
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FIG. 9: (Color online) Elastic energy loss per unit length 
dE/dz of (a) a gluon or (b) a quark in a uniform and static 
medium with a temperature T = 200, 300 and 400 MeV as 
a function of the initial energy Eo from LBT simulations 
of a single scattering (open symbols) or multiple scatter¬ 
ings (L=8fm) (solid symbols) as compared to analytic results 
(solid lines) with a small-angle approximation. 


has a logarithmic dependence on the initial parton en¬ 
ergy. As shown in Fig. [9] (solid lines), this analytic result 
under small-angle scattering approximation agrees well 
with LBT simulations of a single scattering. 

The logarithmic energy dependence of the elastic en¬ 
ergy loss can also explain the difference between single 
and multiple scatterings in Fig. [9] During multiple scat¬ 
terings in the LBT model, energy and momentum is con¬ 
served during each scattering. The leading parton will 


lose its energy along its propagation path. The reduced 
energy of the leading parton will then lead to smaller 
elastic energy loss for subsequent scatterings according 
to the logarithmic energy dependence. This will lead 
to an overall reduction of the averaged energy loss per 
unit distance during multiple scattering over a long dis¬ 
tance. This reduction on the averaged energy loss per 
unit distance is more significant when the total energy 
loss A E = LdE/dz is comparable to parton’s initial en¬ 
ergy. Since a gluon’s energy loss is 9/4 times larger than 
that of a quark, this degradation of energy and the energy 
loss is more significant for a gluon with lower energies in a 
QGP medium at higher temperatures as shown in Fig. [9] 
To illustrate this degradation of energy loss over time, 
we show in Figs. [To] and |TT] the averaged elastic energy 
loss per unit distance as a function of propagation time 
for different initial parton energies and medium temper¬ 
atures. It is clear that the degradation of parton energy 
loss with time is the strongest for a gluon with the lowest 
initial energy and in a medium at the highest tempera¬ 
ture. For an initial high energy quark, the energy loss 
remains almost a constant even at a later time. 


In the LBT model, energy and momentum are con¬ 
served in each parton-medium scattering (with the en¬ 
ergy and momentum of “negative” partons subtracted 
from the sum) along the path of the parton propaga¬ 
tion. The energy lost by the leading parton as shown 
in the above calculations will be carried by thermal re¬ 
coil partons, minus that of “negative” partons, in the 
form of jet-induced medium excitations. To illustrate the 
transfer of energy from the leading partons to jet-induced 
medium excitations, we plot in Fig.[l2]the energy spectra 
of all partons including both leading and the thermal re¬ 
coil partons (with the “negative” partons subtracted) at 
different times (t = 2, 6, and 10 fm/c) during the prop¬ 
agation of an initial gluon with energy E 0 = 100 GeV 
in a uniform QGP medium at a constant temperature 
T = 400 MeV. At the beginning of the parton propaga¬ 
tion, leading partons peak around the initial energy Eq 
with power law tails towards small values due to the dom¬ 
inant i-channel scattering, while thermal recoil partons 
(minus “negative” partons) peak around their thermal 
energy ~ T+fP D /T with a power law tail. At later times, 
leading partons continue to lose energy due to multiple 
scattering and peak at smaller energies ~ E 0 — tdE/dz 
with increasingly broadened distributions. The number 
of jet-induced medium partons (thermal recoil partons 
minus “negative” partons), on the other hand, will con¬ 
tinue to increase with softened spectra which resemble 
a thermal distribution at later times. In principle, the 
spectra of leading and jet-induced medium partons will 
eventually merge and approach a thermal spectrum as a 
result of the equilibration. But this limit is beyond the 
validity of the LBT model. 
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FIG. 10: (Color online) Energy loss per unit path length 
dE/dz of a gluon with different initial energies in a uniform 
medium with a constant temperature (a) T = 200 MeV and 
(b) T = 300 MeV as a function of the propagation time. 


IV. JET-INDUCED MEDIUM EXCITATIONS 

During the parton propagation in a QGP medium, the 
energy and momentum deposited in the medium should 
dissipate with time and propagate like a sound wave. 
Since the leading massless parton travels at the speed 
of light which is much faster than the maximum sound 
velocity c s = 1 /v/3 in an ideal QGP, the sound wave 
induced by the propagating parton becomes supersonic. 
It can form a Mach-cone shock wave as shown in many 
recent studies [i5U55] , Simulations of jet propagation 


FIG. 11: (Color online) The same as Fig. [l0]for a quark. 


within parton transport models also show similar Mach- 
cone-like features [571 Eg. Studies of such jet-induced 
medium excitations can help us to understand the exper¬ 
imental data on dissipation of energy lost by quenched 
jets [71] and extract bulk medium properties such as 
sound velocity or the equation of state (EoS) of the QGP 
medium. Since this article is devoted to a general descrip¬ 
tion and test of the LBT model in a uniform medium, we 
will focus on the effect of dynamical scatterings of prop¬ 
agating partons in a uniform medium. This is similar to 
the effect of viscosity on the propagation of shockwaves 
in a QGP medium [S31Z5HZ7]. We will discuss the prop¬ 
agation of jet-induced Mach cones in an expanding QGP 
medium in the future when we carry out realistic studies 
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FIG. 12: (Color online) Energy spectra of leading partons and 
jet-induced medium partons (thermal recoil partons minus 
“negative” partons) at different times during the propagation 
of an initial gluon with energy Eo = 100 GeV in a uniform 
QGP medium at a constant temperature T = 400 MeV. 

of jet-medium interaction in high-energy heavy-ion colli¬ 
sions. 


A. Diffused Mach cone 

To study jet transport and jet-induced medium exci¬ 
tations at the same time in a QGP medium, both of the 
outgoing partons in each 2 —>• 2 scattering are tracked in 
LBT Monte Carlo simulations. They are allowed to un¬ 
dergo further propagation and scattering in the medium. 
Of the two outgoing partons, the more energetic one is 
designated as the leading parton for further propagation 
and scattering. The flavor of the outgoing leading par- 
ton can be different from incident parton though this is 
quite rare since the scattering is dominated by small- 
angle t-channel processes. The less energetic or soft par- 
ton from each parton-medium scattering is considered as 
the recoil medium parton that is also allowed to further 
propagate and interact with the medium. In addition, 
the initial parton from the thermal medium participat¬ 
ing in the scattering is also tracked and allowed to further 
propagate in the LBT model whose four-momentum is 
subtracted from the final result to account for the back- 
reaction in the Boltzmann transport equation. These 
thermal partons from the back-reaction are denoted as 
“negative” partons in our study and are included in all 
final results on parton spectra and jet reconstruction. In 
general, we refer to the thermal recoil and “negative” 
partons as jet-induced thermal partons. 

Since jet-induced thermal partons are allowed to prop- 


(a) t=4 fm/c 



(b) t=8 fm/c 



FIG. 13: (Color online) Energy density distribution of the jet- 
induced medium excitation by a gluon with an initial energy 
Eo = 100 GeV along (2) and transverse (r) to the initial 
parton direction after (a) 4 and (b) 8 fm of propagation in 
a uniform QGP medium at a constant temperature T = 400 
MeV. The jet parton starts propagating along the 2 direction 
at 2 = 0 and r = 0. 


agate in the QGP medium according to the Boltzmann 
transport equation in the LBT model, we can study the 
linear response of the medium to jet-medium interac¬ 
tion along the path of a jet’s propagation. Interactions 
among jet-induced medium partons are neglected. Such 
an approximation of linear response is valid when the 
jet-induced medium excitation is much smaller than the 
thermal background, <5/ <C /. 

The energy and momentum lost to the medium by the 
propagating parton will act like a source for jet-induced 
medium excitation through the propagation and trans- 
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port of the jet-induced medium partons. Shown in Fig. |13| 
are two-dimensional and contour energy density distribu¬ 
tions of the medium excitation induced by a propagating 
gluon with an initial energy Eg = 100 GeV along (z) and 
transverse (r) to the initial parton direction after 4 fm/c 
(upper panel) and 8 fm/c (lower panel) of propagation 
in a uniform QGP medium at a constant temperature 
T = 400 MeV. The jet parton starts propagating along 
the z direction at z = 0 and r = 0. Note that in each 
event of parton propagation there are only a limited num¬ 
ber of jet-induced medium partons. Figure [j~3| shows the 
energy density distributions of jet-induced medium par- 
tons or medium excitation averaged over many events. 
The number of medium recoil and “negative” partons in¬ 
creases separately with time in LBT Monte Carlo simula¬ 
tions. The net energy density, however, remains approx¬ 
imately constant because of energy-momentum conser¬ 
vation and the near constant energy-momentum transfer 
rate via parton-medium interaction. One can see clearly 
the formation and propagation of a Mach-cone-like shock 
wave trailing the leading parton. The edge of the shock 
wave travels at a speed limited by the velocity of light 
while the peak of the shock wave propagates at the ef¬ 
fective sound velocity in the medium. The shock wave 
is significantly diffused during its propagation because 
of the dissipation due to finite values of viscosities as a 
result of parton-parton collisions as implemented in the 
LBT model. The energy density of jet-induced medium 
partons is negative behind the propagating parton. This 
is the diffusion wake induced by the jet-medium interac¬ 
tion which essentially depletes the thermal medium par- 
ton density behind the path of the propagating parton. 

The diffused shock wave in the shape of an ellip¬ 
tic paraboloid induced by jet-medium interaction within 
the LBT model as shown in Fig. [13] is sharply differ¬ 
ent from the Mach-cone shape given by other linear re¬ 
sponse theories [53], hydrodynamic simulations [5TH55] or 
anti-de Sitter/Conformal field theory (AdS/CFT) corre¬ 
spondence studies isa EH with a source term that trav¬ 
els along a given direction. The leading parton in the 
LBT Monte Carlo simulations, however, does not travel 
along a fixed direction. It receives transverse momentum 
transfers in random directions and accumulates trans¬ 
verse momentum broadening along its propagation path. 
Such momentum broadening amounts to a smearing in 
the direction of the source term for jet-induced medium 
excitation and the width of the smearing increases with 
the propagation time. This eventually leads to the ellip¬ 
tic paraboloid shape of the shock wave instead of a sharp 
Mach-cone. 


B. Azimuthal diffusion 

To illustrate consequences of the transverse momen¬ 
tum broadening of leading partons on the final spectra 
of jet-induced medium partons, we show in Fig. [14] an¬ 
gular distributions of jet-induced medium partons with 
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FIG. 14: (Color online) Angular distributions of medium par- 
tons with energy E = 1.0—1.25 (solid), 1.5-1.75 (dashed), 3.0- 
3.25 GeV (dotted) induced by a gluon with an initial energy 
Eo = 100 GeV after (a) t = 2, (b) 6 and (c) 10 fm/c of propa¬ 
gation in a uniform QGP medium at a constant temperature 
T = 400 MeV. 
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different ranges of energy relative to the initial direction 
of a gluon with energy Eq = 100 GeV after t = 2 (up¬ 
per panel), 6 (middle panel), and 10 fnr/c (lower panel) of 
propagation in a uniform medium at a constant tempera¬ 
ture T = 400 MeV. In a uniform medium, event-averaged 
parton distributions should be azimuthal symmetric with 
respect to the initial direction of the propagating parton. 
For the convenience of future study of azimuthal distri¬ 
butions in a cylindrical frame along the beam direction 
in high-energy heavy-ion collisions, we project the an¬ 
gular distribution onto a plane aligned with the initial 
propagating parton, 


dN 

dcj) 



8 {4> — arctan[tan 9 cos ip ]), 


(18) 


where 6 and ip are the polar and azimuthal angles, re¬ 
spectively, relative to the initial jet parton direction. We 
refer to the angle 4 > as the projected azimuthal angle 
and the distribution as the projected azimuthal distri¬ 
bution. Soft jet-induced medium partons (£1 = 0 — 1 
GeV) are mostly “negative” partons from the the back- 
reaction whose contributions to the spectra are negative 
according to our subtraction scheme. These “negative” 
partons have a broad angular distribution that does not 
change much over time. Contributions to more energetic 
jet-induced medium partons from “negative” partons de¬ 
crease with time, especially along the initial direction of 
the propagating parton, and thermal recoil partons dom¬ 
inate. Intermediate energy thermal recoil partons during 
the early time have an angular distribution similar to 
that of a single scattering as shown in Figs. [5] and [6] that 
has two peaks away from the initial direction. At later 
times, however, the double-peak structure disappears due 
to transverse momentum broadening of the leading par- 
ton and diffusion of thermal recoil partons due to multi¬ 
ple scatterings. The effect of the transverse momentum 
broadening on the angular distribution of thermal recoil 
partons is stronger for less energetic initial partons at 
later times, especially when the total energy loss becomes 
sizable relative to the initial parton energy. 


V. RECONSTRUCTED JETS IN MEDIUM 

In high-energy heavy-ion collision experiments, recon¬ 
structed jets have become a powerful tool for the study 
of jet quenching [T4UT6| 1291432] [751 . These jets are col¬ 
limated clusters of hadrons and are reconstructed from 
the calorimetric energy of final hadrons within a jet cone, 

V(v - Vj) 2 + O - </v) 2 < R, (19) 

using a jet finding algorithm [62], where r) ( r/j ) and (f> 
{4>j) are the pseudorapidity and azimuthal angle of the 
final hadron (jet), respectively. In heavy-ion collisions, 
one should also subtract the background energy within 
the jet cone from hadrons in the underlying events. Such 
background subtraction is extremely nontrivial and may 
affect the final jet energy and jet spectra. 


If we neglect nonperturbative effects of parton 
hadronization, we can approximate final hadronic jets 
with partonic jets reconstructed from final partons in 
theoretical simulations. In our LBT Monte Carlo sim¬ 
ulations, we assume a perfect subtraction of the back¬ 
ground from underlying events if the initial propagating 
parton is not present. We achieve this by including both 
leading partons as well as jet-induced medium partons 
for the jet reconstruction. We use a modified version of 
the FASTJET code [62] with the anti -hr algorithm for 
jet reconstruction in our study within the LBT model 
in which a “negative” parton is treated as a normal one 
when its pseudorapidity and azimuthal angle are calcu¬ 
lated, but its energy and momentum are subtracted from 
the total energy and momentum of the reconstructed jet 
in each iteration of the jet-finding algorithm. We only 
consider the leading jet per event through the FASTJET 
jet finding package with the given algorithm. 


A. Jet energy loss and pr broadening 

During the propagation of a fast parton, it under¬ 
goes multiple scatterings with medium partons in the 
LBT model. The leading parton loses energy and ex¬ 
periences transverse momentum broadening. The lost 
energy and momentum transfer are then carried by jet- 
induced medium partons which will propagate and go 
through further scattering with the medium leading to 
jet-induced medium excitations. Some fraction of the 
energy and momentum carried by jet-induced medium 
partons will still be confined within the jet cone. How¬ 
ever, an increasing fraction will be transported outside 
the jet cone via large-angle scattering and diffusion of 
jet-induced medium partons over the course of the par- 
ton propagation. This will lead to a reduction of the 
energy of the reconstructed jet. Shown in Fig. [15] is the 
energy loss of reconstructed jets with a cone size R = 0.3 
as a function of time for an initial gluon with energy 
Eq = 50 and 100 GeV in a uniform QGP medium with a 
constant temperature T = 400 MeV. The solid symbols 
represent the energy loss of jets reconstructed with all of 
jet-induced medium partons while the open symbols rep¬ 
resent jets without “negative” partons. The jet energy 
loss increases linearly with time or distance of propaga¬ 
tion similarly to the elastic energy loss of a propagating 
parton. 

Th effect of “negative” partons from the diffusion wake 
is negligible during the early stage of the parton propa¬ 
gation because there are very few of them. At the later 
stage, as the number of these “negative” partons becomes 
large, they deplete significantly the thermal medium be¬ 
hind the propagating parton and effectively modify the 
background underlying the jet. When such modification 
of the underlying background is taken into account via 
the subtraction of “negative” partons from the jet cone, 
the effective jet energy loss becomes bigger as shown in 
Fig. [15] Without subtraction of “negative” partons, the 
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FIG. 15: (Color online) Jet energy loss in a uniform QGP 
medium at a temperature T = 400 MeV as a function of 
propagation time of an initial gluon with energy Eq = 50 
and 100 GeV. Jets are reconstructed with all partons (leading 
and jet-induced medium partons) (solid symbols) or without 
“negative” partons (open symbols). 


jet energy loss tends to saturate at later times. The lin¬ 
ear time dependence of the total effective jet energy loss 
is restored only after the energy of “negative” partons is 
subtracted. 

During the propagation of an energetic parton, the 
leading parton experiences transverse momentum broad¬ 
ening through interaction with medium partons. The 
thermal recoil parton in the same scattering should also 
carry the same amount of transverse momentum as the 
leading parton but in the opposite direction. If a re¬ 
constructed jet contains both the leading parton and the 
thermal recoil partons (minus the “negative” partons), its 
total transverse momentum should not be influenced by 
the scattering. However, reconstructed jets do not con¬ 
tain all jet-induced medium partons because of the finite 
jet-cone size R. Large-angle scatterings and diffusion of 
jet-induced medium partons should lead to a reduction 
of the jet energy, as shown in the above, as well as a net 
transverse momentum with respect to the initial parton 
direction. Shown in Fig. [TH] are the transverse momen¬ 
tum distributions of reconstructed jets with respect to 
the initial direction of a propagating gluon with energy 
Eq = 100 GeV in a uniform medium at a constant tem¬ 
perature T = 400 MeV. There is an apparent transverse 
momentum broadening at later times during jet trans¬ 
port, similar to the broadening of leading partons. There 
is also a power law tail in the transverse momentum dis¬ 
tribution. This is caused mainly by a single large angle 
parton-medium scattering during the jet propagation. 
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FIG. 16: (Color online) Jet transverse momentum distribu¬ 
tions at different times during the propagation of an initial 
gluon with energy (a) Eq = 50 GeV and (b) Eq = 100 GeV in 
a uniform QGP medium at a temperature T = 400 MeV. Jets 
are reconstructed with all partons (leading and jet-induced 
medium partons) (solid symbols) or without “negative” par- 
tons (open symbols). 


B. Medium modification of jet structures 


Reconstructed jets in LBT Monte Carlo simulations in 
this study contain both leading partons and jet-induced 
medium partons inside the jet cone. To examine the com¬ 
position of reconstructed jets and their evolution with 
time, we show parton distributions within reconstructed 
jets as functions of the longitudinal momentum fraction 
zj = Pi/Ejet in Fig- Poland z 0 = pl/E 0 in Fig. [l 8 j where 
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Pl is a parton’s longitudinal momentum along the direc¬ 
tion of a reconstructed jet with energy Ej et . Since we 
start with a single energetic parton, the initial parton dis¬ 
tribution function should be a 5 function at Zj = Zq = 1. 
Jet-induced medium partons from parton-medium inter¬ 
action and leading partons with reduced energy populate 
the distribution at zj < 1 or Zq < 1 after the jet transport 
begins. Further scatterings of both leading partons and 
jet-induced medium partons lead to the enhancement of 
soft partons at zj <C 1 or zq <C 1 as shown in Figs. [17] 
and |18| Since the definition of momentum fraction zj is 
normalized by the reconstructed jet energy Ej et which is 
still dominated by the energetic leading parton, the par- 
ton distribution in zj continues to have a peak at zj ~ 1 . 
The parton distribution in zq which is normalized by the 
initial parton energy Eq is, however, suppressed at large 
Zq ~ 1 and its initial peak at large zq ~ 1 disappears at 
later times as the leading parton continues to lose energy 
during its propagation through the medium. Similar be¬ 
havior was reported in the study of y-triggered jets in 
high-energy heavy-ion collisions |60| . The modified frag¬ 
mentation functions in terms of momentum fraction of 
the original jet energy (or the 7 ’s energy) are much more 
sensitive to parton energy loss than that in terms of the 
momentum fraction of the reconstructed jet energy. Since 
“negative” partons are subtracted from both the parton 
distribution and the reconstructed jet energy (reducing 
.Ejet), they will reduce the soft parton distributions at 
small Zj -C 1 or Zq <C 1 and enhance the parton distri¬ 
bution at large zj ~ 1 . 

To study the evolution of the jet transverse structure, 
we show in Fig. [l9]the transverse profile of reconstructed 
jets at different times for a propagating gluon with an 
initial energy Eq = 50 GeV (upper panel) and Eq = 
100 GeV (lower panel) in a uniform QGP medium at a 
constant temperature T = 400 MeV. The jet transverse 
profile is defined as the average fraction of jet energy 
inside an annulus in the 77 — <j> plane, 


p{r) 


1 1 ^ E{r - Ar/ 2 , r + Ar/2) 

A rN& ^ E(0,R) 


( 20 ) 


as a function of r = \J (p — rjj) 2 + ((f> — 4 >j) 2 , where 
E(r\,r- 2 ) is the summed energy of all partons in the an¬ 
nulus between radius r\ and r 2 inside the jet cone. When 
a high energy parton propagates through the plasma, it 
loses energy to the medium. The lost energy is carried 
away by thermal recoil partons whose further diffusion 
through multiple scatterings will transport some of the 
lost energy outside the jet cone leading to the jet en¬ 
ergy loss and broadening of the transverse profile toward 
the edge of the jet cone, as seen in Fig. [19] The “nega¬ 
tive” partons are mostly soft partons at large angles away 
from the initial parton’s direction. When their energy is 
subtracted from the reconstructed jet, the broadening of 
the transverse profile will be reduced, especially near the 
edge of the jet cone. 
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FIG. 17: (Color online) Parton distributions in the momen¬ 
tum fraction zj = Pl/E j e t within the jet cone R = 0.3 at 
different times for a propagating gluon with initial energy (a) 
Eq = 50 and (b) 100 GeV in a uniform QGP medium at a 
temperature T = 400 MeV. Solid symbols represent all par- 
tons including leading, recoil, and “negative” partons while 
the open symbols show results without “negative” partons. 


VI. SUMMARY AND DISCUSSIONS 


In this paper we have described in detail the basic el¬ 
ements of the LBT model for the study of jet transport 
in a QGP medium with the complete set of elastic 2 —> 2 
scattering processes. We verified the Monte Carlo code 
in the LBT model by direct comparisons to semi-analytic 
results of total and differential scattering rates via nu¬ 
merical integrations. We calculated the transverse mo- 
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FIG. 18: (Color online) Parton distributions in the momen¬ 
tum fraction zq = Pl/Eq within the jet cone R = 0.3 at 
different times for a propagating gluon with initial energy (a) 
Eq = 50 and (b) 100 GeV in a uniform QGP medium at a 
temperature T = 400 MeV. Solid symbols represent all par- 
tons including leading, recoil, and “negative” partons while 
the open symbols show results without “negative” partons. 


FIG. 19: (Color online) Jet transverse profile of the leading 
jet by a propagating gluon with initial energy (a) Eq = 50 
and (b) 100 GeV at different times in a static and homo¬ 
geneous medium at temperature T = 400 MeV. Results that 
include all (leading+positive+negative) partons are presented 
in solid symbols while results without “negative” partons are 
presented in open symbols. 


mentum broadening and elastic energy loss of the leading 
parton along the path of the parton propagation. Within 
the pQCD approach to each elastic scattering between 
the leading and thermal medium parton, both the jet 
transport coefficient q and the energy loss per mean-free- 
path have a logarithmic dependence on the parton en¬ 
ergy. Such an energy dependence leads to some nontriv¬ 
ial time or distance dependence of the average transverse 
momentum broadening and energy loss per unit distance 


due to the decrease of the leading partons’ energy over 
time. This dependence is more significant when the total 
energy loss becomes comparable to the initial parton’s 
energy. 

We have also illustrated the dissipation of the en¬ 
ergy and momentum lost by a propagating parton in 
the medium via transport of jet-induced medium partons 
(thermal recoil partons and “negative” partons from the 
back-reaction) in the LBT model. The transport of these 
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jet-induced medium partons within the LBT model effec¬ 
tively forms a supersonic shock wave and a diffusion wake 
behind the leading parton. The supersonic wave has a 
diffused Mach-cone shape because of the transverse mo¬ 
mentum broadening of the leading parton along its path. 
The energy spectra of jet-induced medium partons in the 
supersonic wave are shown to resemble that of a thermal 
distribution at later times of the jet propagation. 

Using a modified version of the FASTJET jet-finding 
package with the anti-/cT algorithm [62| in which the en¬ 
ergy of “negative” partons from the back-reaction is sub¬ 
tracted, we show the effect of jet-induced medium par- 
tons and their dissipation in medium on reconstructed 
jets. Transport of these jet-induced medium partons 
outside the jet cone constitutes a significant reduction 
of jet energy loss which should be taken into account 
for any realistic study of jet suppression in high-energy 
heavy-ion collisions. Their inclusion in the jet reconstruc¬ 
tion also influences the transverse momentum broaden¬ 
ing, parton distributions (fragmentation functions) and 
transverse profiles of reconstructed jets. 

Transverse momentum distributions of both leading 
partons and reconstructed jets are shown to have a power 
law tail underneath Gaussian distributions from multi¬ 
ple scattering. These power law tails come from a sin¬ 
gle large angle parton-medium scattering during the jet 
propagation. Study of these power law tails in heavy-ion 
collisions can shed light on the microscopic nature of the 
QGP medium at different scales. 


For a realistic description of jet transport in the 
QGP medium in high-energy heavy-ion collisions, we 
will have to implement inelastic processes such as gluon 
bremsstrahlung induced by multiple parton-medium 
scattering [sum and couple LBT Monte Carlo simu¬ 
lations with bulk medium evolutions from 3+1D hydro- 
dynamic model calculations, which should be constrained 
by bulk hadron spectra from existing experimental data. 
These parts of the LBT model and the phenomenologi¬ 
cal study of jet suppression in heavy-ion collisions will be 
discussed in subsequent publications. 
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